Weak localization of bulk channels in topological insulator thin film 
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Weak antilocalization (WAL) is expected whenever strong spin-orbit coupling or scattering conies 
into play. Spin-orbit coupling in the bulk states of a topological insulator is very strong, enough 
to result in the topological phase transition. However, the recently observed WAL in topological 
insulators seems to have an ambiguous origin from the bulk states. Starting from the effective 
model for three-dimensional topological insulators, we find that the lowest two-dimensional (2D) 
bulk subbands of a topological insulator thin film can be described by the modified massive Dirac 
model. We derive the magnetoconductivity formula for both the 2D bulk subbands and surface 
bands. Because with relatively large gap, the 2D bulk subbands may lie in the regimes where the 
unitary behavior or even weak localization (WL) is also expected, instead of always WAL. As a 
result, the bulk states may contribute small magnetoconductivity or even compensate the WAL 
from the surface states. Inflection in magnetoconductivity curves may appear when the bulk WL 
channels outnumber the surface WAL channels, providing a signature of the weak localization from 
the bulk states. 

PACS numbers: 03.65.Vf, 73.20.-r, 73.25. -l-i, 85.75.-d 



I. INTRODUCTION 



Topological insulators are materials with gapped bulk 
states but gapless surface states>i"— Due to the topolog- 
ical origin and Dirac fermion nature of the topological 
surface states, the topological insulators are expected to 
have excellent performance in transport 4.i^ One of in- 
triguing transport features of topological insulators is the 
weak antilocalization (WAL), appearing as the negative 
magnetoconductivity with a sharp cusp in low fields i^^— 
WAL is intrinsic to topological insulators: (i) so far, most 
samples have low mobility and long coherence length, 
making the quantum interference an important correc- 
tion to the diffusion transport; (ii) due to the spin- 
momentum locking resulted from strong spin-orbit cou- 
pling, the single gapless Dirac cone of the topological sur- 
face states carries a it Berry phasef2ii2^ which changes the 
interference of time-reversed scattering loops from con- 
structive to destructive. The destructive interference will 
give the conductivity an enhancement, which can be de- 
stroyed by applying a magnetic field that breaks the tt 
Berry phase, leading to the negative magnetoconductiv- 
ity with the cusp. 

WAL is always expected in systems with either strong 
spin-orbit scattering or coupling, -3^'^P Due to vacancies 
and defects, most as-grown topological insulators have 
also bulk band carriers which actually dominate the 
transport i3A'?3 The bulk bands of topological insulator 
possess strong spin-orbit coupling, strong enough to in- 
vert the normal band structure and give rise to the topo- 
logical phase transition that defines the nontrivial na- 
ture of topological insulators. One of the simplest choices 
when considering the conduction bands of a band insu- 
lator with spin-orbit coupling that respects time-reversal 



symmetry is the Rashba model 



Hr = 



2m* 



+ \{(7xky - aykx), 



(LI) 



where h is the Planck constant over 2tt, {kx,ky) is the 
wave vector, k'^ — k^ + ky, m* is the effective electron 
mass, and cr^.y are Pauli matrices that describe the elec- 
tron spin. The Rashba model describes two branches of 
conduction bands, each band has spin locked to momen- 
tum, just like the gapless surface states of topological 
insulator. As a result, the Berry phase for each band of 
the Rashba model also gives exact tt.^^ According to the 
TT Berry phase argument, if the bulk bands of a topolog- 
ical insulator were described by the Rashba model, they 
should also have the weak antilocalization in the quantum 
diffusion transport, just like the gapless surface states. 

Experimentally, the weak antilocalization is studied by 
fitting the magnetoconductivity with the Hikami-Larkin- 
Nagaoka formula,^'' 
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where ^I^ is the digamma function, B is the magnetic field, 
a and C are two fitting parameters, C is an effective phase 
coherence length, a is a prefactor, each band that carries 
a TT Berry phase should give an a = —1/2 prefactor i^ii^ 
In the weak interband coupling limit, multiple indepen- 
dent bands with WAL should add up to give bigger neg- 
ative a, e.g., -1, -1.5. The experimentally fitted a covers 
a wide range between around -0.4 and -1.1, suggesting 
that the observed WAL can be interpreted by considering 
only one or two surface bands,—"— despite of the coexis- 
tence of multiple carrier channels from bulk and surface 
bands at the Fermi surface. Even, the sharp WAL cusp 
can be completely suppressed by doping magnetic impu- 
rities only on the top surface of a topological insulator.— 



Why the observed WAL seems to be weakly tied to the 
bulk bands with strong spin-orbit coupling still poses a 
mystery to both experimentalists and theorists. Various 
interpretations were proposed, such as electron-electron 
interaction^^ and the mixture of the surface states on the 
top and bottom surfaces of the topological thin film.— 

In this work, we try to investigate the role played by 
the bulk states in the quantum diffusion transport of a 
topological insulator thin film. We find that the two- 
dimensional (2D) modified Dirac model can provide a 
unified description for both the surface bands and the 
lowest 2D bulk subbands of a topological insulator thin 
film. We derive the magnetoconductivity formula for the 
2D modified Dirac model in the weak interband scatter- 
ing limit. In this unified description, whether one has 
weak antilocalization or weak localization (WL) is gov- 
erned mainly by the mass (gap) term in the modified 
Dirac model. In the massless limit, one has weak antilo- 
calization, while a finite gap can lead to the weak local- 
ization or the unitary behavior. Contrast to the gapless 
surface states, the bulk states of a topological insulator 
actually have relatively large gap. Therefore, while the 
surface bands probably exhibit weak antilocalization, we 
suggest that bulk bands may reside in the weak local- 
ization or the unitary regime. The experimentally ob- 
served "weak antilocalization" may be a collective result 
from both the weak antilocalization of the surface chan- 
nels and weak localization (or unitary behavior) of the 
bulk channels. This may help to explain why the fit- 
ting parameters of the Hikami-Larkin-Nagaoka formula 
cover a wide range in the experiments.—"— The unitary 
behavior (with small magnetoconductivity proportional 
to the square of the magnetic filed) of bulk channels may 
help to understand why the bulk states seem "missing" 
in the magnetoconductivity although they actually dom- 
inate the transport. 

The paper is organized as follows. In Sec. |TT1 we intro- 
duce the 2D modified Dirac model that describes both 
the surface bands and the lowest 2D bulk subbands of 
a topological insulator thin film. In Sec. IIII[ the mag- 
netoconductivity formula is presented for the 2D mod- 
ified Dirac model in the quantum diffusion regime and 
the weak interband scattering limit. In Sec. IIV Al we 
show the crossover from WAL to WL for a single channel 
of the modified Dirac model. In Sees. IIVBI and IIV CI 
we present the total magnetoconductivity of two surface 
WAL channels and multiple bulk WL channels. Finally, 
a summary is given in Sec. [V] 



II. UNIFIED DESCRIPTION OF BULK AND 
SURFACE OF A TOPOLOGICAL INSULATOR 

The minimal model to describe a three-dimensional 
(3D) topological insulator is the modified Dirac modeli^ 
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FIG. 1: (Color online) (a) The gapped bulk (grey area) and 
gapless surface (dashed lines) bands of a 3D topological in- 
sulator, (b) The quantum confinement along the z direction 
splits the 3D bulk bands into 2D subbands, while the hy- 
bridization of the top and bottom surfaces opens a gap (A) 
for the gapless surface bands, (c) In the ultrathin limit, the 
Fermi surface intersects with only one pair of bulk subbands 
(with band gap M) and one pair of gapped surface bands 
(each curve is two- fold degenerate). The horizontal dot line 
marks the Fermi energy Ef measured from the Dirac point. 



where k — {kx, ky, k^) are wave vectors, the 4x4 Dirac 
matrices ex. = {ax, ay,az), and /3 satisfy the relations 

a^ = j3 — \, UiUj = —oijOii, UiP = —/3ai. 

fc2), M), = m- B{kl + fc2) 



(2.2) 



Sk = C + D(/s2 + fc^), Ml, = m- B{ki + fc^), k± 
kx ± iky, and k"^ — k^ + ky. A, B, C, D, and m are 
model parameters.'^^''^^ For simplicity, three-dimensional 
isotropy is assumed. The modified Dirac model consti- 
tutes the minimal description of the nontrivial topologi- 
cal nature of band insulators: for mB > 0, the solution 
of d-dimensional topologically-protected in-gap bound- 
ary states (surface states, edge states) can be found by 
solving the (d-f l)-dimensional modified Dirac model with 
open boundary conditions.— i^i^i^ Contrast to it, we 
emphasize that the non-Dirac models with spin-orbit cou- 
pling or scattering [e.g., the Rashba model in Eq. (jl.fp ] 
are not enough to describe a topological insulator because 
they can not give the topologically-protected boundary 
states solution. Therefore, it is unconsidered to expect 
the bulk states of topological insulator to have WAL by 
simply assuming them as non-Dirac electron gases with 
strong spin-orbit coupling or scattering. 

With translational symmetry in all three dimensions, 
i?3D in Eq. (|2.II) gives two energy bands [grey area in 
Fig. [la)] 



E: 



3D 
±k 



= ek± JX2^A2(fc2 + A:2), (2.3) 



HsD = Ek + Ak • a + A^k/3, 



(2.1) 



separated by the band gap 2m. Now we consider the film 
geometry by imposing the open boundary conditions to 
the top and bottom surfaces defined along the z direc- 
tion, topologically-protected surface states will emerge in 
the gap [red dashed lines in Fig. [If a)]. For infinite thick- 
ness along the z direction (the bulk limit), the surface 



bands are gapless. Now if we consider finite thickness 
along the z direction, the quantum confinement will split 
the 3D bulk bands E±\^ into a series of 2D subbands 
[solid curves in Fig. [IJb)] as fcz is quantized into dis- 
crete values. Meanwhile, a finite-size gap will open for 
the surface bands [dashed curves in Fig. [Ijb)] due to 
the hybridization of top and bottom surf aces. ^^"''^ In the 
ultrathin limit (e.^., three quintuple layers of Bi2Se3 or 
Bi2Te3 thin film o^^'^'^i"^^'^'^ ), the Fermi surface may inter- 
sect with up to one pair of 2D bulk subbands [see Fig. 
[He)], and the gap A of the surface bands becomes quite 
visible 4^ 



A. Topological surface bands 



properties as fc — >■ oo and give well-defined integer Chern 
number— and Z2 index.— A set of parameters for this 
model is given in Table |l] for 5 nm. 



TABLE I: The parameters in Eq. (|2.5p for d = 5 nm, cal- 
culated from the model parameters of the e ffective model for 
3D topological insulators.— kc = \/M/2B. 





BiaSes 


BiaTeg SbaTes 


D (eV-A^) 


30.4 


49.68 -10.78 


M (eV) 


-0.5 


-0.58 -0.28 


B (eV-A^) 


-44.5 


-57.38 -48.51 


A (eV-A) 


3.33 


2.87 3.40 


kc (A-i) 


0.075 


0.071 0.054 



For the gapped surface bands (dashed curves in Fig. 
[1]), their effective Hamiltonian can be derived from the 
3D model in Eq. (I2T]).3^'^1 

Hs = Dk" + T^{-- Bk^)a, + A{a^ky - ayk,). (2.4) 

It has two 2x2 blocks, with r^ = ±1 the block index. 
'^x,y,z are the Pauli matrices. The model parameters C, 
D, A, B, and A are functions of the thickness. The gap 
A for the surface bands is opened at the Dirac point due 
to the top-bottom surface hybridization. In general, A 
increases with decreasing thickness, but could also van- 
ish at some critical thicknesses where there may be topo- 
logical phase transitions between the quantum spin Hall 
and normal statesi^"— i^ A and B become negligible 
for large enough thickness (e.g., tens of nanometers^). 
A can acquire extra correction in the presence of mag- 
netic dopin g^^i^^ A — >■ A -I- Tz A^, where Aj? represents a 
mean-field from the exchange interaction with magnetic 
impurities. Because block r^ = 4-1 and —1 have opposite 
spin definitions, the same magnetic doping will increase 
the gap in one block while decrease in the other 4i 



In the absence of the Dk'^ and Bk^az terms and for 
only one block of the modified Dirac model, the quan- 
tum diffusion transport is knownJ^ii^ In the gapless 
limit (A/ = 0), there is only weak antilocalization with 
negative magnetoconductivity cusp4S>^ A finite gap M 
will lead to a crossover from weak antilocalization to: 
(i) the unitary regime with small magnetoconductivity 
proportional to the square of magnetic field for mod- 
erate gap/Fermi energy ratios; and to (ii) weak local- 
ization with positive magnetoconductivity cusp for large 
gap/Fermi energy ratiosJ^ii^ 

In contrast to the surface bands with no or small gaps, 
the 2D bulk subbands of a 3D topological insulator have 
relatively large band gaps (more than hundreds of meV, 
see Table. [l|. As a result, instead of always weak an- 
tilocalization, the unitary behavior or weak localization 
is likely expected for the bulk states of a topological in- 
sulators. Moreover, Dk^ and Bk^a^ will bring more sce- 



III. MAGNETOCONDUCTIVITY FOR 2D 
MODIFIED DIRAC MODEL 



B. 2D bulk subbands 

The simplest way to consider the lowest 2D bulk sub- 
bands in Figs. [1] (b) and (c) is to replace (fc^) = and 
{k1) = (tt/c?)^, where d is the thickness of the film. Af- 
ter defining C + Di-n/df = 0, A//2 = m - B{T:/df, 
the Hamiltonian of the lowest 2D bulk subbands can be 
written as 

M 
H = Dk^+ r,(y - Bk^)az + A{a,k, + ayky). (2.5) 

Tz, '^x,y,z have the same meanings as in Eq. p.4[) . Haniil- 
tonians (|2.4I) , (|2.5p , and the BHZ model for the quantum 
spin Hall effect of the HgTe quantum well^ can be clas- 
sified as the modified Dirac model in two dimensionsj^l 
Comparing with the original Dirac model, it has an ex- 
tra Bk^Gz term, which helps to regulate the boundary 



We have shown that the 2D modified Dirac model in 
Eq. (|2.5|) provides a unified description for (i) the 2D 
bulk subbands due to the quantum confinement, and (ii) 
the surface bands of a 3D topological insulator. The 
quantum diffusion transport of the bulk subbands and 
surface bands can be studied starting from this model. 
It has two 2x2 blocks with block index t^ = ±, each 
block has a conduction band and a valence band, de- 
noted as \tz,c/v), with the band dispersions given by 
E^/^ = Dk^ ± ^{M/2-Bk'^Y+A^B, where c and 
V stand for conduction and valence bands, respectively. 
Without loss of generality, we assume that the transport 
is only contributed by the conduction bands that cross 
the Fermi surface. Each band contributes one channel. If 
we assume that the interband coupling by impurity scat- 
tering is weak, the total conductivity is the summation 
of that from each band intersecting the Fermi surface. 
The conduction bands |-l-,c) and | — ,c) should have the 



same magnetoconductivity formula. 

of | + , c) is given by |+, c) = (cos ^ 

where tamp = ky/k^, and at the Fermi surface 



The wavefunction 



2 ,sm^e 



COS0 = 



M/2 



Bkl 



Ef - Dkl ' 



(3.1) 



where Ep and kp are the Fermi energy and Fermi wave 
vector, respectively. 

The weak localization and weak antilocalization hap- 
pen in the quantum diffusion regime, where the mean 
free path {£e) due to elastic scattering is much shorter 
than the sample size, but the phase coherence length {£^) 
due to inelastic scattering is comparable with the sam- 
ple size. The quantum diffusion transport can be studied 
by the standard diagrammatic technique ^^i^iiiSi^S WL 
and WAL are most evident in the magnetoconductivity, 
which is defined as the conductivity change as a function 
of an applied magnetic field B. In this work, we consider 
that the conductivity is measured along the x direction 
and the magnetic field is applied only along z direction. 
In the calculation, we consider the nonmagnetic elastic 
scattering by static centers as well as random scattering 
by magnetic impurities. The calculation is similar to that 
for the topological surface states with a magnetically- 
doped gapfi^ where cos 6 = A/2Ep compared to Eq. 
p.ip . so we present only the result here. For band |-t-, c), 
the zero-temperature magnetoconductivity is found as 



Aa{B) = J2 



4=0,1 



nh 



(1 1 b2 



(3.2) 



where ^ is the digamma function, l\ = fi,/(4e|i3|) is the 
magnetic length, l/^^j = 1/^^ -I- 1/^f , £4, is the phase 
coherence length when cos = 0, 



ao 



2(1 + 1/51)^ 
.90 



"1 



2(l + l/go + l/.g2)' 
51 

2^2(1 + II g^) -i - 2^2(1 + II g^ + II g^) ■ 



n-2 _ yO n-2 

^0 - nl)2{^ , 1 I.. \ *1 



.4^54 ^ip 

go = 2' 



5i=2[ 
92 = 2[ 



a4 (1/^2 + 1/^2) 
1/^2 



(1/^2 _i/^2)(2„252)/(„4+ 54) _ 2/^2 

a^ + b^ lie ^. 



fe4 (1/^2 + 1/^2) 



-1], 

(3.3) 



a = cosf, 6 = sinf. Also, l/£2 = l/£2 + 1/^2^^ ^here 
^e is the mean free path. The calculation assumes short 
mean free path (£e <C ^s). tm is the magnetic scattering 
length that characterizes the strength of the random scat- 
tering by magnetic impurities, l/£^ = l/i1 + l/(.'^-\-l/l1, 
ix,y,z are components along the x, y, z directions, respec- 
tively. In-plane isotropy {^x = (-y) is assumed in Eq. 
(13. 2p . Shorter £„ means stronger magnetic scattering. 

In Eq. dSSl), Tlv=[l- a^b^eiil - e/£l)/{a^ + b'^)] "' 
comes from the correction to velocity from the ladder 



diagrams,^"* and rjn = — (1 — l/r;„ — (? j&^j^ comes from 
the dressed Hikami boxes. ^"^ The lengthes are related to 
characteristic times by ^2 — pj-^ and (^ — Pr^, where 
Tg is the elastic scattering time, t,„ is the magnetic scat- 
tering time. The defined diffusion coefficient T) — v\t 12^ 
where 1/t = 1/re + 1/t„i and the Fermi velocity vp is 
given by 

WF = -ri2DkF 



Asme-2BkpcosQ). (3.4) 



IV. RESULTS 



A. Single channel 
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FIG. 2: (Color online) (a) Magnetoconductivity of a single 
channel of the 2D modified Dirac model. A crossover from 
weak antilocalization to weak localization is expected when 
increasing cos0. £^ = 300 nm and £e ~ 10 nm. Weak 
isotropic magnetic scattering is assumed, £x/V3 — £z/V^ ~ 
im ~ 10000\/3 nm. (b)-(g) Possible pseudospin patterns of a 
conduction band of the modified Dirac model. 

Wc have shown that for each channel contributed by 
the modified Dirac model, the magnetoconductivity is 
given by Eq. p.2p . Despite of many parameters in Eq. 
p.3p , the most decisive one is cos Q defined in Eq. p.ip . 
As shown in Fig. [2J the magnetoconductivity formula 
describes a crossover from weak antilocalization to weak 
localization, controlled by cos0. In the limit cos 8 = 0, 
one has only WAL with negative magnetoconductivity 
cusp. Increasing cosO will drive the system first into 
the unitary regime with small magnetoconductivity pro- 
portional to i?2, and finally to WL with positive mag- 
netoconductivity cusp when cos0 — >■ 1. In the limit 
cos 8 = 0, one has ag = and ai = —1/2; In the limit 



COS 8 -^ 1, one has ao = 1/2 and ai = 0, corresponding 
to that the prefactor a in the Hikami-Larkin-Nagaoka 
formula goes from —1/2 to 1/2. 

The TT Berry phase explained the weak antilocaliza- 
tion of a single gapless Dirac cone,^^ the crossover as 
the function of cos 9 can be understood similarly. For 
the conduction band |+, c), the Berry phase at the Fermi 
surface is found as 



Two surface channels and two bulk channels 



7r(l-cose). (4.1) 



2^ Q 



where cosO is given by Eq. (|3.1|) . Therefore, one has 
the TT Berry phase for WAL when cos 6 = and Berry 
phase for WL when cos 6^1. Different from the Berry 



phase in absence of Bk az term,'^ one may also have 
the TT Berry phase as well as the weak antilocalization at 
2ttp = kp = M/2B, where p is the sheet carrier density 
of investigated bandJ^ Moreover, Dkp may effectively 
reduce Ep to Ep — Dkp, reminding that the effective 
model for the gapless surface states also contains the Dk^ 
term,^^'^^ this may explain the earlier crossover to weak 
localization at relatively large Epi^ Because the Berry 
phase depends on cos 6, the crossover can also be under- 
stood by the pseudospin orientation along z direction. 
Several typical orientation cases are given in Figs. [2Ib)- 
(g). WAL with a = —1/2 corresponds to purely in-plane 
polarization, as in Figs. [2jc) and (e); WL with a = 1/2 
corresponds to the fully polarized along z, as in Figs. [21 
(b) and (g). Due to the Bkpaz term, one may expect 
purely in-plane polarization at kp = y^M/2B, and the 
pseudospin polarization can change from z to ^z direc- 
tion as shown in Figs. [2Kd)-(f). Due to the Dkp, one 
may expect the full polarization along z at finite kp as 
shown in Fig. [DJg), instead of always at the bottom of 
the band (Fig. ^). 

The Berry phase of the valence bands is given by 
7r(l -t- cos 9). Therefore, if the sample is p-type, i.e., the 
Fermi surface is at the 2D valence subbands, the quantum 
diffusion transport will give the similar competition be- 
tween WL and WAL, as long as the 2D valence subbands 
are also well described by the massive Dirac model (sit- 
uations in Bi2Se3 and Bi2Te3 may be more complicated 
because of their valence band maxima away from the F 
point). 

Recently the quantum diffusion transport was also the- 
oretically studied for the HgTe quantum wellfi^ where 
the suppression of weak antilocalization (corresponding 
to the unitary behavior in this work) was also found for 
the BHZ model. Beyond the suppression of the weak an- 
tilocalization, we also expect the weak localization in the 
HgTe quantum well because the BHZ model can also be 
classified as the modified Dirac model. 

Because the two blocks of the modified Dirac model 
in Eq. (|2.5p are mutually time-reversal partners, the 
crossover here does not break time-reversal symme- 
try, different from those in the ferromagnetic semicon- 
ductors with spin-orbit couplin g^^i^^ and magnetically- 
doped topological surface states i^'^'^^ 
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FIG. 3: (a) The magnetoconductivity of two surface chan- 
nels with cos 0s = (dotted line), two bulk channels with 
cos Qb ~ 0.85 (dashed line), and their summation (solid line), 
(b) solid: the same as the solid line in (a). Dash line: fit- 
ting to solid line between B € [—0.5, 0.5] with a — —0.316 
and jC — 941. Dotted line: fitting to solid line between 
B G [-0.05, 0.05] with a = -0.659 and C = 363. All the 
parameters are given in Fig. (2] 
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FIG. 4: Total magnetoconductivity of two surface channels 
in the weak antilocalization regime and two bulk channels in 
the unitary or weak localization regime. All the parameters 
are given in Fig. [21 

Now we consider the case shown in Fig. [ijc), with 
two channels of surface bands and two channels of bulk 
subbands (note that each curve in Fig. [T] represents two 
degenerate bands labeled by the block index t^). We 
have shown that for each channel of the modified Dirac 
model, the magnetoconductivity is given by Eq. p.2p . 
which describes a crossover from weak antilocalization 



TABLE II: The Fermi energies of n-type as-grown bulk crys- 
tals and thin films of topological insulators from the ARPES 
measurements. The Fermi energies are measured from the 
bottom of the conduction band. 



Ref. 


sample 


thickness Ferm 


I energies (eV) 




BizSeg 


bulk 


0.1 


31 


BiaSea 


bulk 


0.05 


^ 


BiaTeg 


bulk 


0.045 


is 


BiaSea 


bulk 


0.15 


42 


BiaSea 


2-5 QL 


0.1~0.2 


43 


BiaTea 


80 nm 
2-5 QL 


0.03 
0.1-0.2 


13 


BiaSea 


SQL 


0.11 



TABLE III: Fitting the weak locahzation curves in Fig. |4]by 
the Hikami-Larkin-Nagaoka formula in Eq. (|1.2p . a and L 
are two fitting parameters. B are the perpendicular magnetic 
field in units of Tesla. 





BG 


[-0.5, 0.5] 


Be 


-0.05, 0.05] 


cos 0s 


cosOfl 


a 


L (nm) 


Of 


£ (nm) 





0.5 


-0.883 


341 


-0.991 


301 





0.7 


-0.650 


468 


-0.944 


309 





0.8 


-0.444 


683 


-0.815 


331 





0.85 


-0.316 


941 


-0.659 


363 





0.9 


-0.184 


1502 


-0.419 


433 


0.05 


0.5 


-0.787 


193 


-0.945 


169 


0.05 


0.7 


-0.489 


266 


-0.830 


175 


0.05 


0.7 


-0.119 


582 


-0.315 


221 


0.1 


0.5 


-0.594 


117 


-0.829 


98.5 


0.1 


0.7 


-0.211 


181 


-0.505 


108 



to weak localization with increasing cos0. In the weak 
interband scattering limit, the total magnetoconductivity 
is the summation of that from each band. Similar to Eq. 
TTj) . we introduce two control parameters 



cos9r = 



cos 6 c = 



A'f/2 - Bk\ 
Ef - Dkj, ' 
A/2 - Bkj. 



Ep 



Dkl 



(4.2) 



hand, we expect relatively large cos 9b that gives the 
unitary behavior or weak localization. When the Fermi 
surface is near the bottom of the band, the gap M and 
the Fermi energy Ep dominate cosQb- The bulk sub- 
bands have relatively large gap, about 0.5 eV for Bi2Se3 
and Bi2Te3 (note that this gap is measured at the F 
point, '^^ not the smaller indirect gaps which are irrelevant 
to the topological natur o^^i^^'^^ ). The Fermi surfaces of 
as-grown samples are usually less than 0.20 eV measured 
from the bottom of the bulk conduction band (see Ta- 
ble |ll|, i.e., Ep is about 0.25 to 0.5 eV measured from 
the Dirac point. According to these data, cos &b ranges 
between 0.5 and 1, i.e., between the unitary and weak 
localization regimes according to Fig. [2l 

In Fig. ^a.), we show the total magnetoconductivity 
for two surface and two bulk channels, with cos 65 = 
and cos Qb — 0.85. Although the surface and bulk chan- 
nels exhibit the weak antilocalization and the weak lo- 
calization, respectively, they collectively behave like the 
weak antilocalization. Therefore, in Fig. ^h), we try to 
fit the total magnetoconductivity by the Hikami-Larkin- 
Nagaoka formula. Although widely exploited in fitting 
experiments, this formula may not be suitable for multi- 
ple channels, because the magnetoconductivity consists 
of two competing terms even for single channel of the 
modified Dirac model, not to mention multiple chan- 
nels. Therefore, no ubiquitous fitting result of a and 
C can be reached when changing the fitting range. For 
B e [-0.5,0.5], we obtain a = -0.316 and £ = 941, and 
for B e [-0.05,0.05], a = -0.659 and C = 363. Both 
fittings can not reproduce the magnetoconductivity well 
in the whole range. 

In Fig. m we show more general cases with cos 65 £ 
{0,0.05,0.1} and cos 65 G {0.5,0.7,0.8,0.85,0.9}. They 
are fitted by the Hikami-Larkin-Nagaoka formula in Eq. 
(|1.2p . Note that the fitting results almost cover the ex- 
perimentally recorded a (g [-1.1,-0.38], see Table ITV| . 
These fitting results imply that one can approach to 
a — —1/2 anyway, because of the freedom in choosing 
C and fitting range. Therefore, a ~ —1/2 does not al- 
ways mean that there is only one surface channel of single 
massless Dirac cone. It may happen to be the collective 
result of multiple surface and bulk channels. In Figs. 
SJb) and (c) , one even has weak localization, which is ex- 
pected by us in some ultrathin films though so far there 
was no such report. 



to characterize the bulk subbands and the surface bands, 
respectively. The parameters in Eq. (|4.2p were defined 
in Hamiltonians (|2.4|) and (j2.5p . Note that the Fermi en- 
ergies Ep and Ep are measured from the Dirac points, 
which may be defined differently for the bulk and surface 
bands. In absence of magnetic impurities, two surface 
channels have the same cos 65, and two bulk channels 
have the same cos 63. Except in the ultra thin limit 
(usually less than 5 nm in thickness), the surface states 
has ignorable A and jj jiii38,42^43 j-j- jg jr^^jj. ^^ expect small 
cos 8s that gives weak antilocalization. On the other 



C. Two surface channels and more than two bulk 
channels 

Now we considered the case in Fig. [lib), where more 
than two bulk channels can participate in the transport, 
and the surface channels have negligible gap. Without 
loss of generality, we consider four bulk channels. It is fair 
to expect the 2D bulk subbands higher than the lowest 
pair to have similar descriptions as the modified Dirac 
model in Eq. (I2.5p but with different parameters. For 



TABLE IV: The experimentally fitted prefactor a and phase 
coherence length L in the Hikami-Larkin-Nagaoka magneto- 
conductivity Eq. (|1.2p . 



Ref. 


sample 


thickness 


r(K) 


a 


L (nm) 


^ 


BizSes 


10 nm 


L8 


-0.5~-0.38 


106-237 " 


^ 


BiaSes 


10~20 nm 


0.3 


-0.38 




io 


BiaTea 


50 nm 


2 


-0.39 


331 


11 


BiaSes 


2~6 QL 


1.5 


-0.39 


75~200 


12 6 


BiaSes 
(Bi,Pb)Se3 


45 nm 


0.5 


-0.31 
-0.35 


1100 
640 


ii 


BiaSea 


5~20 nm 


0.01-2 


-l.l~-0.4 


143~ oo 


15 


BiaSea 


3~100 QL 


1.5 


-0.63~-0.13 


100~1000 


i6 


BiaSea 


20 nm 


0.3~100 


-1.1—0.7 


15~300 



"Calculated from B^ = ft/{4e£"'). 

''Electron-electron interaction was also considered in the fitting. 
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0.4- 


"" 




^ 
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FIG. 5: (a) The magnetoconductivity of two surface channels 
with cos 0s = (dotted line) and four bulk channels with 
cos 0s = 0.85 (dashed line), and their summation (solid line), 
(b) Total magnetoconductivity of two surface channels with 
cos 0s = and four bulk channels in the unitary or weak 
localization regime. All the other parameters are given in 
Fig.H 



simplicity, we assume that the four bulk channels have 
the same cos 0b ■ 

Because now we have unequal numbers of weak an- 
tilocalization channels from the surface bands and weak 
localization channels from the bulk subbands, the weak 
localization behavior of the bulk channels may be hidden 
by the weak antilocalization of the surface channels at 
small magnetic field, but the superiority in channel num- 
ber of the bulk band will eventually change the trend at 
some finite magnetic field. As shown in Fig. [51 inflections 
in magnetoconductivity appear at finite magnetic fields, 



and even change the magnetoconductivity from negative 
to positive. These infiections will be evident signatures 
indicating that the bulk states could have weak localiza- 
tion instead of always weak antilocalization. 



V. SUMMARY AND DISCUSSION 

Starting with the modified Dirac model as a unified 
starting point, we derive the magnetoconductivity for- 
mula for both the lowest 2D bulk subbands and surface 
bands of a topological insulator thin film. A crossover 
from the weak antilocalization to weak localization is ex- 
pected, controlled by the pseudospin polarization defined 
by the parameters in the modified Dirac model. Unlike 
those in the ferromagnetic semiconductor, this crossover 
does not break time-reversal symmetry. Due to their rela- 
tively large gap, we suggest that the bulk states may lie in 
the regimes where the unitary behavior or even the weak 
localization is also expected, instead of the expected weak 
antilocalization. As a result, the experimentally observed 
weak antilocalization may be a collective result of the 
weak antilocalization of the surface bands and weak lo- 
calization of the bulk bands. It may explain the deviation 
of the fitting prefactor of the Hikami-Larkin-Nagaoka for- 
mula from the expected —1/2 in the experiments, as well 
as the insensitive response of the bulk states to magnetic 
field or doping. When the bulk channels outnumber the 
surface channels, infiection in the magnetoconductivity 
may appear at some finite magnetic field. The inflection 
will give a signature that the bulk states of a topological 
insulator can give weak localization, although they have 
strong spin-orbit coupling. 

In this work, we consider only the ultrathin limit of 
a 3D topological insulator. In the bulk limit, because 
a large number of 2D bulk subbands will contribute to 
the transport, the direct coupling [by the k^ terms in 
Hamiltonian (12.11) ] and scattering via impurities among 
them become inevitable. In this limit, we argue that 
the inter-subband coupling and scattering will give ex- 
tra "Cooperon gaps" to all the vertex corrections from 
the maximally crossed diagrams for the 2D bulk bands. 
These Cooperon gaps will reduce eflfectively the phase 
coherence length, and drive the transport via each 2D 
bulk subband out of the quantum diffusion regime. As a 
result, we expect the 3D bulk states to have neither weak 
localization nor weak antilocalization, but the unitary be- 
havior with small magnetoconductivity proportional to 

We thank H. T. He, J. N. Wang, F. C. Zhang, Junren 
Shi, W. Q. Chen, and B. Zhou for helpful discussions. 
This work is supported by the Research Grant Council 
of Hong Kong under Grant Nos. HKU 7051/lOP and 
HKUST3/CRF/09. 
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